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ABSTRACT 

In this paper, we explore how the forthcoming generation of large-scale radio contin- 
uum surveys, with the inclusion of some degree of redshift information, can constrain 
cosmological parameters. By cross-matching these radio surveys with shallow optical 
to near-infrared surveys, we can essentially separate the source distribution into a low- 
and a high-rcdshift sample, thus providing a constraint on the evolution of cosmolog- 
ical parameters such as those related to dark energy. We examine two radio surveys, 
the Evolutionary Map of the Universe (EMU) and the Wcsterbork Observations of the 
Deep APERTIF Northern sky (WO DAN). A crucial advantage is their combined po- 
tential to provide a deep, full-sky survey. The surveys used for the cross-identifications 
are SkyMapper and SDSS, for the southern and northern skies, respectively. We con- 
centrate on the galaxy clustering angular power spectrum as our benchmark observ- 
able, and find that the possibility of including such low redshift information yields 
major improvements in the determination of cosmological parameters. With this ap- 
proach, and provided a good knowledge of the galaxy bias evolution, we are able to 
put strict constraints on the dark energy parameters, i.e. wq = —0.9 ± 0.041 and 
w a = —0.24 ± 0.13, with type la supernovas and CMB priors (with a one-parameter 
bias in this case); this corresponds to a Figure of Merit (FoM) > 600, which is twice 
better than what is obtained by using only the cross-identified sources and greater 
than four time better than the case without any redshift information at all. 

Key words: large-scale structure of the universe - cosmological parameters - cos- 
mology: observations - radio continuum: galaxies. 



1 INTRODUCTION 

A new era for precision cosmology with radio continuum 
surveys is close to becoming a reality. Several experiments 
across the world are currently beginning operations, such as 
the LOw Frequency ARrajQ (LOFAR) and APERTIF0 at 
the Westerbork Synthesis Radio Telescope (WSRT), in the 



* E-mail: stefano.camera@ist.utl.pt (SC); mgrsantos@ist.utl.pt 
(MGS); david.bacon@port.ac.uk (DJB). 

1 http://www.lofar.org 

2 http:/ /www. astron.nl/general/apertif/apertif 



northern hemisphere, and the projects for the Australian^ 
and African]] SKA Pathfinders, in the southern hemisphere. 

Recently, iRaccanelli et all (|2012t) have shown promis- 
ing forecasts for constraining cosmological parameters using 
some of the large-scale radio continuum surveys planned for 
these forthcoming telescopes. In this paper, we investigate 
the usefulness of including redshift information, for the sub- 
sample of radio continuum sources which are successfully 
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cross-identified with optical surveys containing photomet- 
ric or spectroscopic redshifts. We examine the impact of 
this partial redshift knowledge on cosmological constraints 
with forthcoming radio surveys. The underlying effect of this 
multi-wavelength information is reasonably simple to un- 
derstand: shallow optical surveys should be able to identify 
most of the low redshift galaxies from these large radio sur- 
veys, allowing us to "extract" the high redshift tail expected 
for the radio distribution of sources. 

Although the detailed analysis is more involved, basi- 
cally the combination of these two "macro" redshift bins 
will give us a powerful constraint on the redshift evolution 
of cosmological models as an approximate alternative to de- 
manding redshift surveys. For the purpose of this work and 
in order to clarify the improvements achieved with this red- 
shift information, we will focus on the radio galaxy two- 
point correlation function as the cosmological observable; 
there are several other available probes such as the inte- 
grated Sa chs- Wolfe effect and cos mic magnification — as dis- 
cussed in iRaccanelli et all (|2012f ) — which will merit further 
investigation in the context of specifi c modified matter and 
gravity theories (|Camera et alj|2012l ). 

The paper is structured as follows. We briefly present 
the EMU and WODAN surveys and their main characteris- 
tics in § O while the surveys used for cross-identifications, 
providing redshift information, are described in §[3] §|4] out- 
lines the data analysis procedure which we will consider for 
obtaining cosmological constraints. In §[5] we give the most 
important formulae for the galaxy clustering angular power 
spectrum and its statistics. In § [U we introduce the fidu- 
cial dark energy (DE) model used in our analyses and in 
§ [7] we describe our method. In § [8] we present the analysis 
of the combined radio and optical surveys for constraining 
the DE equation of state. First, we present the cosmological 
constraints available when there is no redshift information 
at all; secondly, we show what would be obtained in the 
idealised case of full knowledge of every radio-source red- 
shift. Finally, we describe the realistic scenario expected to 
be available in the recent future, where we cross-identify as 
many EMU/ WO DAN galaxies as possible with optical red- 
shift surveys, in particular SkyMapper and SDSS. For this 
last case, we scrutinise a simple one-parameter model for the 
radio-galaxy bias, as well as a two-parameter model which 
accounts for a further redshift dependence. Conclusions are 
drawn in §[§] 



2 FORTHCOMING RADIO SURVEYS 

Here, we present the two survey designs we use in our cal- 
culations. A summary of their properties is given in Ta- 
ble [1] where N g represents the total number of detected 
radio galaxies at IOct and z m the median redshift (see 
IRaccanelli et al.ll2012l for details). It is worth noticing that 
these two surveys have the same redshift distribution of 
sources, shown in Fig. [T] (black histogram). Furthermore, 
they will cover the entire sky, if their data are combined. 
This is one of the major strengths of the present analysis, 
providing a homogeneous all-sky catalogue, since the two 
surveys have basically the same sensitivity, galaxy number 
density and redshift distribution. 



Table 1. Specifics of the EMU and WODAN surveys, where N g 
is the total number of detected sources at 10<7, and z m the median 
redshift. 

Survey Area Frequency Sensitivity N g z m 

EMU 3tt 1400MHz lO^Jy 2.2 ■ 10 7 1.1 

WODAN 7T 1400 MHz 10/fJy 7.3 ■ 10 6 1.1 




z 



Figure 1. Redshift distribution of sources dN/dz{z) for the 
EMU/WODAN surveys. The area is normalised to the mean 
galaxy number density. The red histogram shows the distribution 
of radio objects which have optical counterparts in SkyMapper 
and SDSS. 

2.1 EMU 

EMU (Evolutionary Map of the Universe. iNorris et al1l201ir i 
is an all-sky continuum survey planned for the new Aus- 
tralian SKA Pathfinder CASKAP. lJohnston et aTlbOOSt ) tele- 
scope under construction on the Australian candidate SKA 
site in Western Australia. The primary goal of EMU is to 
make a deep 10 /^Jy /beam radio continuum survey of the 
entire Southern Sky, extending as far North as +30° at a 
resolution of ~ 10 arcsec; it will also have sensitivity to ex- 
tended structures. The EMU survey is expected to begin in 
2013. 



2.2 WODAN 

WODAN (Westerbork Observations of the Deep Apertif 
Northern sky survey) is planned to chart the entire north- 
ern sky above +30° down to a proposed rms flux den- 
sity at 1.4 GHz of 10 ^tjy per be am with a resolution of 
~ 15 arcsec (|Rottgering et al]|201 It will be able to do this 
because of the new phased array feeds (APERTIF) being 
put on the Westerbor k Synthesis Radio Telescope (WSRT, 
lOosterloo et alj l2010h . The current schedule for the com- 
mencement of APERTIF surveys is 2013. 

LOFAR is another experiment which merits comment. 
This is a multi-national telescope with stations spanning Eu- 
rope; the core of LOFAR is situated in the north-east of the 
Netherlands, with stations on longer baselines both within 
the Netherlands and across Germany, UK, France and Swe- 
den. Other stations may also be added throughout the rest 
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of Europe in the coming years. LOFAR large-area contin- 
uum surveys probe to similar depth and source densities 
as EMU and WODAN rtRottgerinel 120031 : iRottgering et all 
l201ll : lvan Haarlem fc et alj|2012h ." and can therefore be used 
on their own or in conjunction with EMU and WODAN for 
cosmological studies. This is because the lower frequency of 
LOFAR makes it sensitive to different source populations 
and provides spectral information. However, for the rest of 
this paper we concentrate on only the EMU and WODAN 
surveys, noting that LOFAR could also be used in the north- 
ern hemisphere. 



3 SURVEYS FOR CROSS-IDENTIFICATIONS 

In this section, we present the surveys we shall use to cross- 
identify the EMU/WODAN sources. Our goal is to ob- 
tain any sort of redshift information for some of the radio 
sources that will be detected by the EMU and WODAN 
surveys using surveys available around the same time. We 
shall mainly focus on current and forthcoming large-scale 
surveys which have multi-filter data across the optical and 
near-infrared to facilitate reasonably accurate photomet- 
ric r edshifts, eithe r through pure templat e fitting methods 
(e.g. lBenitedl2OO0l ; iBolzonella et al.ll200ol ). template fitting 
methods combined w i th known spectroscopic redshifts (e.g . 
iFeldmann et al.ll200rj ; lllbert et aLll2006l ; ISalvato et ai1l201ll ) 
or wi th empirical meth o ds, such as artificia l neur al networks 
(e.g. iFirth et all |2003l: ICollister fc Lahavl |2004|) or Gaus- 



Table 2. Adopted bins for cross-identifications with SkyMappcr 
and SDSS (equal size bins for each range). The assumed scatter 
for photo-2 bins is 0.05(1 + z) for photo-z. 



sian Processes (e.g. Way fc Srivast ava 2006; Bonficl d et alj 
l2010h . 

The most important surveys to co nsider for these ra - 
dio surveys are therefore SkyMa pper (|Keller et al.l 120071 ') 
and SDSS (|Eisenstein et al.1 l201lf) at optical wavelengths, 
combined with the UKIRT I nfrared Deep Sky Surv ey Large 
Area Survey (UKIDSS-LAS; lLawrence et al.ll2007n and the 
proposed UKIRT Hemisphere Survey at northern latitudes, 
along the VISTA Hemisphere Survey at Southern Latitudes. 
We also note that t he Wide-field Infrared Survey Explorer 
l |Wright et alj|20ld . WISE; ) could also play a role in help- 
ing with photometric redshifts over both hemispheres. Fur- 
ther in the future, there are deeper planned surveys ob- 
serving large areas, such as th ose using the Lar ge Synop- 
tic S urvey Telescope (LS ST; llvezic et al] 120081 ) and Eu- 
clicjfl (|Laureiis et al.|[201ll ). We do not consider the slightly 
smaller area survey s such as the Dark Energy Survey (DES; 
lAbbott et aLll2005l ) in this paper, but note that they could 
provide important additional information by constraining 
redshifts to fainter magnitudes and higher redshifts. 

In this paper, we focus on SkyMapper and SDSS. The 
former is a 20, 000 deg 2 optical survey of the southern sky, 
whilst the latter will cover 12, 000 deg 2 in the northern sky. 
For both surveys, we adopt a conservative magnitude limit 
of r < 22 mag which should ensure that the redshifts can 
be measured to reasonable accuracy and certainly to the 
level required for our investigations here. SkyMapper should 
provide photometric redshifts whilst SDSS should be able to 
also provide spectroscopic redshifts — at least for low z's. 

To assess the number of radio sources in both WODAN 
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SkyMappcr 



SDSS 



Spectroscopic 
Photometric 



5 bins (0.0 < z < 0.1 



3 bins (0.0 < z < 0.3) 
3 bins (0.3 < z < 0.8) 



and EMU that will have an optical counte r part t o r ~ 22, we 
use the recent results of iMcAlpine et all (|2012l ). who anal- 
ysed how many reliable cross-identifications to radio sources 
can be found as a function of resolution and depth of the 
radio data coupled with the depth of the optical and near- 
infrared data, using the V ISTA Deep Extrag alactic Obser- 
vations (VIDEO) Survey jjarvis et alj|20ia i. In Fig. Q] we 
show the expected redshift distribution of radio sources de- 
tected at > Wa in EMU and WODAN (black) along with 
the expected number of optical counterparts at r < 22 from 
SDSS and SkyMapper (red). 



4 STRATEGY FOR REDSHIFT BINNING 

In this section, we describe a realistic scenario for an analysis 
leading to cosmological constraints including redshift infor- 
mation. First, we note that the redshift distribution shown in 
Fig. Q] should be known with reasonable accuracy for forth- 
coming radio surveys. This is because there are radio surveys 
such as the Austra lia Telescope Large Area Survey(ATLAS; 
iNorris et ai]|2006D and radio surveys cov ering the UKIDSS- 
Ultra Deep Survey jsimpson et alJl2006T ) and the COSMOS 
field (jSchinnerer et alJl201Cf ). covering a few square degrees 
to the same depth as EMU and WODAN, and for which deep 
spectroscopic data is becoming available. We shall consider 
an analysis where a fraction of the EMU/WODAN galaxies 
are cross-identified with SkyMapper and/or SDSS. Hence, 
for those sources we shall be able to perform an actual red- 
shift binning. 

We assume that spectroscopic redshifts have a negligi- 
ble uncertainty, and that they will be available for the radio 
sources using the SDSS up to z ~ 0.3. For those sources 
without spectroscopic redshifts, namely radio sources with 
identifications in SkyMapper and those without spectro- 
scopic redshifts in SDSS (i.e. with z > 0.3) we shall rely 
on photometric redshifts and assume a photometric-redshift 
scatter a? h = 0.05(1 + 2). The binning we use is summarised 
in Table [2] Note that this does not refer to all the sources 
that will be detected by these optical surveys, but only to 
those that are cross-matched to EMU and WODAN. In par- 
ticular, radio AGNs detected in the radio occupy massive 
galaxies, which should correspond to the SDSS LRGs (Lu- 
minous Red Galaxies), so in principle we could have spec- 
troscopic information from SDSS to higher redshifts (e.g. up 
to z ~ 0.5 instead of z < 0.3 that we assume here). 

On the other hand, one of the most promising properties 
of EMU and WODAN is the high-redshift tail of their source 
distribution (see Fig. [TJ, which is composed of galaxies diffi- 
cult to cross-identify with SkyMapper or SDSS due to their 
faint nature. Nevertheless, these can still bring important 
additional information to the analysis. We will assess the im- 
pact of these objects by making calculations for both a cross- 
identified catalogue, and a catalogue whose source distribu- 
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tion is given by the difference between the EMU/ WO DAN 
source distribution and that of the survey used for the cross- 
identifications, fn addition, we note that the sky coverage 
of the optical surveys will not overlap exactly with EMU 
or WODAN. Therefore, we can also make calculations for 
a third catalogue which follows the EMU/ WODAN source 
distribution and covers all the sky not surveyed by the op- 
tical surveys. 

To summarise, the analysis for this realistic-binning 
strategy will combine three "effective surveys," whose 
dN/dz(z) are defined as follows: 

i) all the EMU/ WODAN galaxies cross-identified by 
SkyMapper/SDSS, appropriately binned in redshift (see 
Table 0; 
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ii) all the EMU/ WO DAN galaxies which are in the 
same patch of sky as SkyMapper/SDSS but which are not 
cross-identified ; 

Hi) the part of the EMU/WODAN survey probing the sky 
left uncovered by SkyMapper/SDSS. 



5 GALAXY ANGULAR POWER SPECTRUM 

We now introduce the main cosmological observable consid- 
ered in this work. For our radio surveys, let us consider the 
2-point angular correlation function £ 9 of radio sources. It 
can be seen as the excess probability S g of finding two radio 
sources at a certain physi cal separation s = x — y one from 
each other (|Peebleslll980h , viz. 



(1) 



here, Sd is Dirac's delta function. By performing the Fourier 
transform of such an observable we obtain the radio source 
over-density power spectrum P 9 (k, z), where k is the physi- 
cal wavenumber, with k = |k|. This is related to the under- 
lying total-matter power spectrum P s (k,z) of the matter 
density contrast S = Sp/p by the radio-source bias function 
b g , which we discuss in greater detail in § 15.11 Thus, we have 
P 3 (k,z) = b g 2 P s (k,z). 

However, when we have no radial information, we in- 
stead deal with projected quantities. Thus, the angular 
power spectrum C 3 (£) of the radio source fluctuation — 
where £ is the angular wavenumber — then reads 



C 3 {£) = 4tt [W 3 (£,k)] 2 P 9 (k,z = 0), 



with W 3 (£, k) a proper line-of-site weight function. A widely 
used sim plification is given by the so-called Limber's approx- 
imation ((Kaiser 1992), where £ = k\\ here, x{ z ) ls the radial 
comoving distance, which is related to the Hubble expan- 
sion rate H(z) via dx = dz/H(z). Limber's approximation 
is valid when ^ > 1, but it has been s hown that the conver- 
gence is already good for £ > 10 (e.g. [Hullioooh . Therefore, 
it is a suitable approximation, since for larger angular scales 
the cosmic variance uncertainty is dominant. In this limit 
we have 



C 3 {£) = Jd X 



W 3 ( X ) 



(3) 



Figure 2. The bias multiplied by the redshift distribution of the 
radio sources as a function of redshift for the EMU and WODAN 
surveys (solid line). Also shown are the product of the bias and 
redshift distribution for surveys at 1 mJy, 5 mJy and 10 mJy to 
illustrate the advance that is expected from the new generation 
of radio surveys compared to surveys such as the NVSS. For the 
fu nctional forms of the b ias of each sour ce population, we re fer 
to iRaccanelli et al ] (|2012l . Fig. 3: see also IWilman et~aDl2008l l . 



with W 3 (x) defined through 



dN 



W 3 [ X {z)] = H{z)b g {z) — {z) 



dz 



(4) 



5.1 Radio-Galaxy Bias 



The radio-source bias is worth a deeper digression. In thi s 
paper, we adopt the approa ch of IRaccanelli et al.l (|2012l b 
They used the simu lations of lWilman et all l|2008h (sec also 
IWilman et"aHl20ld ) for the SKA continuum survey. By gen- 
erating source catalogues with the S-cube databastjf] corre- 
sponding to the radio flux-density limits of the proposed 
EMU and WODAN surveys, they obtained an estimate of 
the bias for each source population used here. No t e that 
the source distribution and bias of IWilman et al.l (|2008T ) 
are compatible wi th several observati onal results such as 
the NVSS survey (|Condon et alj|l998h . In Fig. H we show 
the product of the bias and the redshift distribution of the 
sources, i.e. W g (z)/H(z); we also illustrate the products of 
the bias and redshift distribution for surveys at 1 mJy, 5 mJy 
and lOmJy, to illustrate the advance that is expected from 
the new generation of radio surveys compared to surveys 
such as the NVSS. 

In this paper, we assume a reasonable knowledge a pri- 
ori of the bias evolution, so that only a one- and a two- 
parameter model are used, as explained below. The analysis 
of small sky patches with deep radio data and detailed red- 
shift information from other multi-wavelength data should 
provide a reasonable measurement of the bias directly from 
the redshift space distortions without assumptions on the 
DE model, and thus allow for this more simple parameter- 
isation (note that the bias is assumed scale independent). 



http://s-cubcd.physics.ox.ac.uk 
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Nevertheless, as constraints on the DE evolution become 
more stringent, a more precise procedure will be required, 
such as a joint fit to both the bias and DE model using the 
deep, small survey and the large, shallow radio survey com- 
bined. Specifically, from correlation function of the faint ra- 
dio sources from the deep survey, we shall have important in- 
formation on the galaxy clustering as a function of redshift, 
whilst the data from shallow, large-scale surveys will provide 
information on the scale-dependent effects of DE. Moreover, 
this could include other estimators such as the ISW and, in 
particular, the cosmic magnification, which is less sensitive 
to the bias issue. Finally, other possible paramctcrisations of 
the bias within the theoret ical framework of the halo model 
(e.g. ICoorav fc Shethll2002l , for a review) should be explored; 
they could allow us to reduce the number of parameters. The 
current paper addresses the constraining capabilities of the 
next generation of radio surveys that will soon be available, 
when combined with other multi-wavelength data using the 
clustering estimator as a benchmark, and we hope that the 
results obtained open the window for further studies. 

The effective bias b g (z) of Eq. @ accounts for the dif- 
ferent biases for each source population. The main results of 
this work (§ 18. 3[) have been obtained by allowing for the 
amplitude of this effective bias to take any value, while 
maintaining the redshift evolution of lRaccanelli et all |2012l . 
their Fig. 3). That is to say, we parameterise the effective 
bias as A b b g (z), with a fiducial value A b = 1, and we then 
marginalise the whole Fisher matrix over such amplitude. 

Nonetheless, we also study the effects of allowing for a 
greater uncertainty in the redshift evolution of the bias by in- 
troducing another parameter that accounts for an additional 
redshift dependence. Specifically, we make the substitution 



A b bg(z) -»• (A b + B b z) b g {z), 



(5) 



with B b — as a fiducial value. We shall see that this will not 
qualitatively change our results — though the presence of one 
more free parameter necessarily broadens the constraints. 



0.3 
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Figure 3. Relative difference between DE and ACDM galaxy 
clustering angular power spectra. The black (solid) curve refers 
to the unbinncd EMU+WODAN case, the red (short-dashed) 
and the green (long-dashed) curves correspond — according to Ta- 
ble^ — to the first SDSS and the last SkyMapper bin, respectively, 
and the blue (dot-dashed) line corresponds to the EMU/WODAN 
sources without cross- identifications (see Fig.[TJ. 



where a — 1/(1 + z) is the scale factor. Thus, /de(z) in 
Eq. © is 



/de(z) = exp 



-3 I a da' 1 + W ^ 



(8) 



Hereafter, we shall use t he well-known Chevallier-Polarski- 
Lind er parameterisation (|Chevallier fc Polarskill200ll ; lLinderl 
12003? ) 



6 DYNAMICAL DARK ENERGY 

Here, we briefly outline the cosmological model tested in 
our analysis. Since the objective of this work is to test the 
effects of redshift information, we choose a dynamical DE 
model as our benchmark. We consider a flat Friedmann- 
Lemaitre-Robertson- Walker Universe, filled with a perfect 
fluid of baryons, cold dark matter and DE. Their abundances 
in units of the critical density are Sl b , SIdm and Que = 
1 — Q m , respectively, and the total matter contribution is 
Q. m = Q. b + Qom- Therefore, the expansion history of the 
Universe H = dlna/dt is (see e.g. lBartelmannll2010l ) 

H{z) = ^^.(I + ^' + SIde/deW, (6) 

with Ho — 100 /ikm/s/Mpc the Hubble constant and /de(z) 
the (unknown) function governing the time evolution of DE. 

For an exhaustive and recent review on DE in gen- 
er al, we suggest to the read e r the comprehensive monograph 
of lAmendola fc Tsuiikawal (|2O10T l . For the purpose of this 
work, it is sufficient to say that, in principle, there is a time- 
dependent equation of state linking the DE pressure to its 
energy density, namely 

Pde(o) = wde(o)pde(«), (7) 



u)de((i) = w + (1 - a)w a , (9) 

which is a Taylor expansion around a cosmological constant 
(today). Indeed, the ACDM limit is reached when wo — > — 1 
and w a —¥ 0. 

To highlight why including some redshift information 
is important for the aim of constraining DE with the ra- 
dio source angular power spectrum, in Fig. [3] we show the 
C 3 (£) of Eq. © for different redshift bins. Ultimately, this 
means that the angular power spectra can in principle probe 
departures from standard ACDM, if the signal-to-noise is 
sufficiently high. We plot the relative difference between 
the DE and ACDM galaxy-clustering angular power spectra 
as a function of the angular scale £ for the different effec- 
tive surveys introduced in § U Specifically, the solid (black) 
curve refers to the standard, unbinned EMU+WODAN case, 
while the short-dashed (red) and the long-dashed (green) 
curves respectively correspond to the first SDSS and the 
last SkyMapper bin (see Fig.[T]and Table[2]); the dot-dashed 
(blue) line corresponds to the high-z tail of the effective sur- 
vey (ii). In the next section we shall calculate the degree to 
which we can measure these departures from ACDM given 
realistic source densities. 
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7 METHOD AND FISHER MATRIX 
FORMALISM 

The expected errors on the model paramete rs can be es- 
timated via the Fisher information mat rix (|Fisherl 1 19351 ; 
lJungman et~aill996l ; iTegmark et "alii 19971) . This has the ad- 
vantage that different observational strategies can be anal- 
ysed and this can be very valuable for experimental design. 
The Fisher matrix gives the best errors to expect, and is ac- 
curate if the likelihood surface near the peak is adequately 
approximated by a multivariate Gaussian. Hence, the Fisher 
matrix is simply the sum of the Fisher matrices of each £ 
mode, 



a/3 



E 



2£ + l 8C a {l) rrfl ,,.-i dC%l) 
2 /sky M a 1 1 ' d#g ' 



(10) 



where / s k y is the fraction of the sky covered by the survey 
under analysis and i? = {fi m , fib, wo, w a , h, n s , erg, At} is 
the parameter set, with At the overall galaxy-bias amplitude 
(as introduced in § [5TTJ . Cf' 3 is the covariance for a given I 
mode, 



0(4 + K 



(ii) 



where the second term is related to the Poisson noise, N g 
being the mean number density of sources per steradian. 

We remind the reader that we are considering Limber's 
approximation, the sum in Eq. (|10|) therefore starts at £ m i n = 
10. To check the robustness of the results obtained under 
such assumption, we also performed the analysis with £ m i n = 
2. However, we found that the enhancements thus yielded 
are only of a few percent. This is because, at such large 
scales, the cosmic variance heavily spoils the signal. He nce, 
we fin d it reasonable to neglect £ < 10. Furthermore, iHul 
|2000l ) have shown that Limber's approximation leads to 
an overestimation of the signal for I < 10, which means 
that the actual improvement in the analysis would be even 
smaller than what we have estimated. On the other hand, 
imax = 1000 is kept fixed. 

At this stage, the overall galaxy-bias amplitude Ah rep- 
resents a "nuisance" parameter, since this quantity is not 
one of the fundamental cosmological parameters that we 
want to constrain. This results in a general broadening of 
all the marginal errors, particularly of the power-spectrum 
normalisation erg, which is completely degenerate with the 
bias amplitu de. Therefor e, we marginalise over both these 
parameters (|Wang| |2008j) • To avoid any possible numeri- 
cal instability in the marginalising procedure, we calcu- 
late the the Fisher m atrix marginalised over erg and Ab as 
jAlbrecht et al.ll2009l ) 



F^UA _1 U T F^, 



(12) 



where we define (p = {fi m , fib, Wo, w a , h, n s } and ip = 
{erg, At,}; therefore, is the block of the total Fisher ma- 
trix containing the parameters we want to constrain, whilst 
F** is the nuisance-parameter Fisher sub-matrix. Here, A 
is the diagonal matrix whose elements are the eigenvalues of 
F**, whilst U is the orthogonal matrix diagonalising F**. 
By using Eq. (|12p , our marginalising procedure is more sta- 
ble, since degeneracies in F vv are properly propagated to 
G with no instabilities, and we do not even worry about a 



possibly ill-conditioned F** sub-matrix, since we check its 
stability on the fly by the diagonalisation. 

Another important quantity that we calculate is the cor- 
relation between the parameter pair (pt-tpj , which is defined 
afl 

r w ,-^ = /,J =. (13) 



This quantity tells us whether the two parameters are com- 
pletely uncorrelated (if it is zero) or completely degenerate 
(if it is ±1), or have some intermediate level of correlation. 

Finally, we calculate the DE Figure of Merit (FoM), as 
intro duced by the Dark Energy Task Force l|Albrecht et al.l 
120091) 



FoM 



dot (G" 



(14) 



this is proportional to the inverse of the area encompassed 
by the ellipse representing the 68% confidence level in 
the (wq, ui a )-plane. Therefore, the tighter the constraints, 
the larger the FoM. We also combine all our results with 
the CMB Fisher matrix o btained for the Planck survey 
(|Planck collaboration! 120051 ) and the latest measurements 
of type la supernovae (SNela) from the Union2 sample 
i|Amanullah et al.|[201(il ). 

Since one of our main goals is to show the capability 
of our method of discriminating between the fiducial DE 
model and ACDM, there is also another possible approach. 
There is a pivot redshift z p where the uncertainty on the 
equation-of-state parameter is minimised for a given data 
model. Hence, rather than using Wo, we can introduce 



w p = wo — 



1 + z p 



(15) 



It is also easy to demo nstrate that FoM = (<r Wp a Wa ) 1 holds 
|Albrecht et alj|2009l) . 



8 CONSTRAINTS ON DARK ENERGY 

In order to assess the impact of using redshift informa- 
tion in our analysis, we start by considering a flat Uni- 
verse in the matter-dominated era with matter fraction 
fi m = 0.28, baryonic contribution fit = 0.045, Hubble con- 
stant (in units of 100 Mpc -1 ) h — 0.7, slope of the pri- 
mordial power spectrum n 3 — 0.96 and rms mass fluctua- 
tions on the scale of 8 h~ x Mpc erg = 0.8. The fiducial model 
we use is slightly away from the ACDM point in the DE 
parameter space: {w p, w a } = {—0.9, —0.24}, as fitted by 
IZhao fc Zhangl (|2010l ) against a large number of available 
data. The transfer function for the scale dependence of mat- 
ter densi ty perturbations is calc ulated via the fitting for- 
mulae of lEisenstein fc Hul (|l998l ). which take the baryonic 
features in the matter power spectrum into account. 

We use the redshift distribution of sources of the EMU 
and WO DAN surveys (Fig. [TJ), as outlined in S I2.H2.2I The 
median redshift of the survey is z m = LI- The total number 



7 Note that indices a,/3 = 1 ... 8 refer to the complete set of 
parameters whilst i, j = 1 . . .6 label the sub-set <p. 
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of sources is 2.2 • 10 7 , for EMU, and 7.3 ■ 10 6 , for WODAN. 
The total sky coverage is 4ir, three quarters coming from 
EMU and one quarter from WODAN. As pointed out above, 
we should have access to this overall redshift distribution 
beforehand. 



8.1 No Redshift Information 

First, we present the constraints that can be put on the cos- 
mological parameters if no redshift information is available. 
The marginal errors a ipi on the fiducial values (i lpi of the 
parameters tpi, as well as the correlations between the pa- 
rameters and wq or w a , are shown in Table [3] The numbers 
in brackets refer to results including the full Fisher matrices 
for Planck and Union2. These marginal errors are approx- 
imately of the same order of magnitude as the parameter 
fiducial value. However, when we add CMB and SNela infor- 
mation, the forecast significantly improves. Nevertheless, it 
is still impossible to discriminate between the ACDM model 
and DE. This is because at the time of recombination the DE 
component was massively subdominant compared to matter 
and radiation, and the use of CMB is therefore not enough. 
Similarly, though they are a powerful probe for the late-time 
accelerated expansion of the cosmos, SNela are only sensi- 
tive to the background dynamics. 

A possible solution to this problem is to include in 
the analysis infor mation from other cosm ological probes, as 
has been done bv lRaccanelli et al.l (|2012h . Alternatively, we 
can investigate how much information is hidden in the red- 
shift distribution of the radio sources detected by EMU and 
WODAN. 



8.2 Full Redshift Information 

Before showing the results obtained with the realistic strat- 
egy outlined in § [4l we briefly discuss an extreme idealistic 
case for comparison. Here, we assume we perfectly know 
the redshifts of all the EMU and WODAN sources. Then, 
if we used the same bins of width 0.1 as shown in Fig. [1] 
we would have at our disposal 33 spectra combined into a 
tomographic matrix C^(£). Each element of the matrix is 
a galaxy angular power spectrum as in Eq. ([3]), with the 
fraction of sources contained in the corresponding bin. It is 
worth noting that such a tomographic matrix is diagonal, 
since any cross-correlation between different bins vanishes. 
In other words, sources in the i-th bin do not overlap in 
redshift with sources in the j-th bin. 

The marginal errors a ipi on the fiducial values \i lf>i of 
the parameters tpi, as well as the correlations between the 
parameters and wo or w a , are shown in Table 2] The num- 
bers in brackets refer to the results including the full Fisher 
matrices for Planck and Union2. 

In this extreme case, the results are unsurprisingly ex- 
cellent. The possibility of binning the EMU and WODAN 
sources allows us to track the evolution of DE up to high red- 
shifts, thanks to the long tail of the distribution of the radio 
sources (see Fig. Indeed, even though DE is a subdomi- 
nant species compared to matter until very late times, it is 
crucial to study its evolution, in order to constrain w a . More- 
over, the combination of both surveys yields a full-sky cover- 
age, thus providing better statistics. As a result, the forecast 



marginal errors on the DE parameters are a WQ = 0.022 and 
<J-w a = 0.058, which become a wo — 0.012 and a Wa = 0.033 
when combined with Planck and Union2. These constraints 
finally give FoMs of 1.9 • 10 3 and 5.4 • 10 3 , respectively. 

The correlation coefficients r Vi ~ WQ and r lpi - Wa are par- 
ticularly useful to understand whether this present binning 
method is able of disentangling the parameter dependencies 
and lift some of their degeneracies. Indeed, if we look at 
the correlation coefficients in Table |3j we find that some of 
them are near unity in absolute value. This means that these 
constraints are strongly degenerate with each other; this is 
because the closer to unity is the value of r lpi - ipj , the more 
correlated are the two parameters and ipj . On the other 
hand, the correlations obtained in combination with Planck 
and Union2 are far better for all the parameters. 



8.3 Realistic Binning Strategy 

Now, we present the realistic results that the procedure dis- 
cussed in §[4] will yield. Since we are now dealing with the 
cross-identified sources, their redshifts (and the relative error 
we have to assume) are dependent on the optical survey used 
for the cross-identification, as summarised in Table [2] The 
EMU radio sources will be identified by SkyMapper, which 
will provide photometric estimates for the source redshifts. 
Therefore, we divide the distribution dN/dz(z) of Fig. [T] 
(red curve) into five bins up to z = 0.8 and assume a scatter 
of h = 0.05(1 + z). For WODAN, the situation is slightly 
more complicated. This is because SDSS will provide spec- 
troscopic measurements for objects with z < 0.3, whilst we 
will have to rely on photometric redshifts for the sources at 
higher redshift. Hence, in this case we subdivide the source 
distribution into three spectroscopic redshift bins of width 
0.1 up to z — 0.3, plus three photometric redshift bins up 
to z ~ 0.8. As described in § [4] we then consider two more 
bins, but that do not correspond to a specific redshift: one 
for all the radio galaxies in the same sky area covered by 
SkyMapper/SDSS which are not cross-identified (case u of 
the "effective surveys" ) and another for all other radio galax- 
ies in patches of the sky not covered by the redshift surveys 
(case Hi). 

Fig. [4] presents the 1-a confidence region in the 
(wq, w a )-plane as obtained by using the radio sources (black 
ellipse) and when we also add Planck and Union2 (red el- 
lipse); for comparison, we also show the ellipse for Planck 
and Union2 alone (grey ellipse). The marginal errors o lPi 
on the fiducial values y ltpi of the parameters ipi, as well 
as the correlations between the parameters and wo or w a , 
are shown in Table [5] Again, the numbers in brackets refer 
to results including the full Fisher matrices for Planck and 
Union2. 

As expected, these results are not as tightly constrain- 
ing as those obtained in the idealistic case of a perfect 
knowledge of all the source redshifts, but are a substantial 
improvement over the no-redshift case. Indeed, we obtain 
0"iu o = 0.063 and a Wa = 0.195, which become a wo — 0.041 
and a Wa = 0.13 when combined with Planck and Union2. 
The correlation coefficients again help to obtain a deeper 
understanding of these results. As in the idealistic case, the 
tightness of the constraints on the parameters is fortified 
by the modest degeneracies yielded by our method. For this 
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Table 3. Cosmological constraints when assuming no rcdshift information for the EMU/WO DAN sources. Marginal errors a Vi on the 
fiducial values [i Vi of the DE model parameters <pi and their correlation with wo and Ulfl, V"iz. T(p£ — IV q cllld Tip^ — uj a 

; erg and b g have been 

marginalised over. Numbers in brackets refer to the inclusion of Planck and Union2 Fisher matrices. 
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Figure 4. Marginal l-cr contour in the (wq, to a )-plane of the 
DE model for our realistic binning strategy. The black ellipse 
refers to the Fisher matrix obtained by using only EMU and 
WODAN galaxies and their cross-identifications with SkyMappcr 
and SDSS, whilst the red ellipse includes Planck and Union2 pri- 
ors; for comparison, the grey ellipse shows the contour for Planck 
and Union2 alone. 



purpose, the inclusion of CMB and SNela priors is often 
crucial. 

For instance, in Fig. \S\ we show, for some of the cases 
scrutinised, the power spectra C g (£) (solid, black) with the 
corresponding error bars given by 



AC 9 (I) = 



(2£+l)/ sky 



(16) 



in comparison with a DE model whose {wo, w a } values are 
3cr away according to the errors quoted in Table [5] (dashed, 
red). Specifically, clock- wise starting from the top- left panel: 
the power spectrum for the EMU+WODAN sources without 



redshift information; that for the EMU+WODAN sources 
not cross-identified by either SkyMapper or SDSS; the C 9 (l) 
of the fifth, and last, photometric bin of the sources cross- 
identified by SkyMapper; and that of the first spectroscopic 
bin of the radio galaxies cross-identified by SDSS. 

We can quantify how much our approach permits a dis- 
crimination between DE and ACDM by determining the 
FoMs. In our realistic case, the constraints on wo and w a 
yield FoM = 382 or 225 for the radio sources with or without 
CMB and SNela, respectively. All these major results, to- 
gether with the marginal errors on w a and on the maximally- 
constrained value of w)de(z) — w p of eq. (|15[1 — are presented 
in Table g] 

Finally, we want to understand how the assumptions 
made upon the redshift evolution of the bias influence our 
results. To this aim (see § 15. ip we make use of the ad- 
ditional nuisance parameter Bb, as in Eq. (0, thus intro- 
ducing a further redshift dependence of the bias. We again 
perform the same Fisher analysis as before, where the set 
of nuisance parameter we want to marginalise over is now 
■tp — {as, Ab, Bb}. As expected, the presence of one more 
free parameter leads to weaker constraining power. As a re- 
sult, the Wo-Wa ellipse is broader than in the realistic case. 
Nevertheless, we want to demonstrate that, even in this pes- 
simistic case, our approach gives interesting results. 

Therefore, to conclude, Fig. [fj] summarises the primary 
results of this work. It depicts the FoM versus the bin- 
ning strategy adopted, thus showing the impact of our pro- 
posed pipeline. The labels of the horizontal axis refer to 
the unbinned case (unbinned), the case with cross-identified 
sources alone (X-IDs) and the case resulting from the combi- 
nation of all the effective surveys (TOT) — as ex- 
plained in § [4] Black diamonds refer to the use of only EMU 
and WODAN galaxies and their cross-identifications with 
SkyMapper and SDSS, whilst red circles point to the results 
obtained by the inclusion of the Fisher matrices for Planck 
and Union2. Solid lines refer to the main result for the real- 
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Table 5. The same as Table [3] for the realistic binning strategy. 
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Figure 5. Comparison between the angular power spectra for the fiducial DE model (solid, black with errorbars) and a DE cosmology 
with {wq, Wa} values 3tr's away (dashed, red). Top-right panel: EMU+WODAN sources without any redshift information. Top-left 
panel: EMU+WODAN sources not cross-identified by SkyMapper or SDSS. Bottom-left panel: first spectro-2: bin of the WODAN sources 
cross-identified by SDSS. Bottom-right panel: fifth (and last) photo- z bin of the EMU sources cross-identified by SkyMapper. 



istic case, as highlighted in the previous section; dashed lines 
instead show what happens in the case just described of a 
two-parameter bias. It is clear that the inclusion of cross- 
identifications and, furthermore, the high-z radio tail still 
yields a significant enhancement. 



9 SUMMARY AND CONCLUSIONS 

In this paper, we have explored the impact that redshift in- 
formation could have on the cosmological potential of the 
forthcoming generation of large-scale, radio continuum sur- 
veys. Specifically, we have investigated how the possibil- 
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Table 6. Summary of the errors on wo, w p and w a , and the FoMs for the unbinned, idealistic and realistic cases without (or with) 
Planck and Union2. 





No binning 


Idealistic case 


Realistic case 




2.47 (0.11) 


0.022 


(0.012) 


0.063 (0.041) 


cr Wp 


0.29 (0.018) 


0.0090 


(0.0057) 


0.020 (0.013) 




9.60 (0.40) 


0.058 


(0.033) 


0.195 (0.13) 


FoM 


0.40 (140) 


1909 


(5378) 


260 (603) 



6h 



Radio + SNela + CMB 




O Radio only . 










— - 


^^-- — '^ 




unbinned X — IDs 


TOT 


Binning strategy 





Figure 6. FoMs with or without the inclusion of Planck and 
Union2 (red circles or black diamonds, respectively) versus the 
binning strategy adopted: "unbinned" pertains to the case with 
no rcdshift information, "X-IDs" to the EMU/WODAN sources 
cross-identified by SkyMapper/SDSS and "TOT" to the inclusion 
of all the effective surveys. Solid lines refer to the main case, whilst 
dashed lines to the case with an additional linear paramctcrisation 
for the bias. 

ity of binning the redshift distribution of the radio sources 
would improve the constraining power of these surveys. For 
this purpose, as a fiducial model we adopted a dynami- 
cal DE model whose fiducial values have been fitted by 
IZhao fc Zhangl (|2010D as {w , w a } = {-0.9, -0.24}. We 
have focused on two forthcoming wide-field radio surveys: 
EMU and WODAN. Thanks to their sensitivity and the 
specifics outlined in Table [T] they could be combined thus 
yielding a full-sky survey. This is extremely important, since 
it leads to larger sample size and allows us to lessen the im- 
pact of cosmic variance. 

We have performed a Fisher matrix analysis to forecast 
how the combination of EMU and WODAN will be able 
to constrain the cosmological parameters of the DE model, 
particularly focusing on its extra parameters, wq and w a , 
and using the angular power spectrum of radio sources as 
the cosmological probe. We have examined the impact of 
including redshifts for sources cross-identified with optical 
redshift surveys. We have described results when there is 
no knowledge of redshifts, summarised in Table in this 
case the DE equation of state is poorly determined. We 
have also presented the most idealistic case where all the 
EMU and WODAN redshifts are known, in which case we 
have unsurprisingly obtained extremely tight constraints on 
all the model parameters, with marginal errors on the DE 



equation-of-state parameters a wo = 0.012 and cr Wa — 0.033, 
when combined with Planck and Union2. 

Next, we moved to a more realistic scenario. We have 
included forthcoming optical surveys as cross-identifiers 
for the EMU/WODAN sources; specifically, we have used 
SkyMapper for the southern sky, while SDSS provides red- 
shifts in the north. It is clear that not all of the radio galaxies 
that will be detected by EMU and WODAN will be cross- 
identified by other surveys, as these radio surveys cover the 
whole sky and reach very high redshift. To take this into 
account, we have developed a model of the future prac- 
tical analysis, using three "effective surveys" as described 
in § HI For these sub-surveys, the redshift distribution of 
sources will be those of: i) all the EMU/WODAN galax- 
ies cross-identified by SkyMapper/SDSS, suitably redshift 
binned; ii) all the EMU/WODAN galaxies which are in the 
same patch of sky as SkyMapper/SDSS but which are not 
cross-matched; and Hi) an EMU/WODAN survey probing 
only the sky left uncovered by SkyMapper/SDSS. With this 
approach we have obtained a FoM of 603 and 260 for the 
radio sources with and without CMB and SNela priors, re- 
spectively. These results are very competitive when com- 
pared with other surveys that will also be available in the 
near future. The improvement obtained by this cosmolog- 
ical analysis with radio surveys stems from the fact that 
by cross-correlating with shallow "redshift surveys" we can 
make use of the long high redshift tail of the non-identified 
radio galaxies to provide an extra handle on the evolution 
of DE. This effect can be clearly seen in Fig. [1J] where the 
jump from the second to the third point is only due to the 
extra information from the radio-surveys. 

We want to emphasise that we have accounted for the 
source bias amplitude At, by using a nuisance parameter 
which was allowed to vary freely. However, reality might 
be more complicated with effects on the bias such as scale 
dependence and a change in the relative amplitudes of the 
bias for each finer redshift bin. We leave a deeper analysis 
for future work. It is worth noting that the radio source bias 
could be constrained as a function of redshift by deep cross- 
correlation studies with optical data over smaller areas, but 
any remaining uncertainty requiring more bias parameters 
will worsen the FoMs. To simply test whether our approach 
is robust even in the presence of a more complex bias, we 
have also studied a slightly more complex scenario, in which 
we have added a linear correction to the bias amplitude. 
This means that we deal with one more nuisance parameter. 
Nevertheless, Fig.|6]shows that the main results of this work 
still hold, though the constraints get broader, as expected. 
On the other hand, if we parameterised the functional form 
of the bias not only by its amplitude but also taking its shape 
into account, we should be able to lift the strong degeneracy 
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present here between Af, and the matter power spectrum 
normalisation ag. 

We have also not accounted for possible contaminants 
to the power spectrum that could affect our analysis, since 
the objective of this paper is to concentrate on the impact of 
redshift information for cosmological applications. In partic- 
ular, double sources could affect the angular power spectrum 
at smaller angular scales (< 0°.l). requ iring some cleaning 
technique (see e.g. lOverzier et all 12001 for the case with 
NVSS data). Note, however, that for the 10a-10/xJy flux 
cut we are considering here, these double sources should be 
much less abundant and much less dominant than for the 
NVSS case, so that the contamination might be less severe 
on these sub-degree scales; moreover, we are cutting off our 
analysis at ^ max = 1000. 

Besides, instrumen tal effects may also alter the cl uster- 
ing signal. For instance, iBlake. Mauch fc Sadler! (|2004T ) have 
pointed out that, for angular separations around ~ 0°.2, 
there is an unexpected deficit in NVSS data, which is very 
likely an instrumental effect due to the shortest b aselin e used 
in the observations. Similarly, IBlake fc Walll (|2002T ) have 
shown the risk in not properly considering and modelling 
gradients in the surface density. Indeed, such effects may sig- 
nificantly influence the imprint of the large-scale structure. 
This is because, if we tried to fit the data with a theoretical 
angular power spectrum which does not account for these 
systematics, we would reconstruct the systematic density 
gradients as well as the fluctuations due to clustering. How- 
ever , such problems should arise w i th flux cuts > 10 mjy (see 
also Blake. Mauch fc Sad ler 2004; Bla ke. Ferreira fc Borrull 



I2004F T thus leaving our analysis free of this problem. Nev- 
ertheless, we plan to return to the specific issue of contam- 
ination and cleaning of the cosmological signal in a future 
analysis. 

As a final remark, we note that the optical surveys we 
have used for the cross-identifications should be available at 
approximately the same time as EMU and WO DAN, around 
2015. Therefore, we can see that future large-scale radio sur- 
veys can not only provide a different observational wave- 
length window for constraining the cosmological model, but 
also provide competitive DE constraints. Related to this, 
we aim to exploit such a successful method in a forthcoming 
pape r to also test modified gravity models (|Camera et al.l 
I2012T I. 
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